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Schrodinger's equation with the attractive potential V(r) — —Z/(r q + fi q ) 5 , Z > 0, j3 > 0, q > 1, is 
shown, for general values of the parameters Z and (3, to be reducible to the confluent Heun equation 
in the case q — 1, and to the generalized Heun equation in case q — 2. In a formulation with correct 
asymptotics, the eigenstates are specified a priori up to an unknown factor. In certain special cases 
this factor becomes a polynomial. The Asymptotic Iteration Method is used either to find the 
polynomial factor and the associated eigenvalue explicitly, or to construct accurate approximations 
■ for them. Detail solutions for both cases are provided. 
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I. INTRODUCTION 
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Analytic wave functions corresponding to solutions of Schrodinger's time-independent equation H ^ = for standard 
quantum-mechanical systems, such as the hydrogen atom and the harmonic oscillator, have usually been constructed 
by transforming the original equations in to a differential equation whose solutions are known in terms of certain special 
l— ~~ '■ functions. Indeed the general solutions of the resulting Kummer equation corresponding to such model systems, as 
I \ expressed in the form of confluent hypergeometric functions, have led to a deeper understanding of these systems, 
under both free [J 0] and spatially confined conditions 3] . In this paper we consider the model quantum system 
defined by the Hamiltonian (in atomic units m = h = e = 1) 

3 : 1 z 

S: H=--A + V q (r), V q (r) = - — r . (1) 

2 (r« + 01) 1 

c4 ■ 

— 1 ■ The potential V q (r) represents a family of soft-core (truncated) Coulomb potentials, which are useful as model PO- 
O'S I tentials in atomic and molecular physics. The bound states are obtained in terms of three potential parameters: the 
coupling Z > 0, the cut-off parameter > 0, and the power parameter q > 1. The cases q — 1 and 2 are of special 
physical significance [1, [H, H, 0, H, Gil [O, EH ■ The potential V\ represents the potential due to a smeared charge 
and is useful in describing mesonic atoms. The potential V2 is similar to the shape of the potential due to a finite 
nucleus and experienced by the muon in a muonic atom. Extensive applications of the soft-core Coulomb potential, 
V2, have been made through model calculations corresponding to the interaction of intense laser fields with atoms 
O Q EE ES El El • The parameter can be related to the strength of the laser field, with the range — 20 — 40 
covering the experimental laser field strengths [l3| . 

In an earlier paper (l9j , we have carried out a general analysis of the characteristic features of the energies and wave 
functions of the complete family of soft Coulomb potentials defined by V q . The Schrodinger operator H is bounded 
below. This may be shown immediately by an application of the operator inequality (20l . [2l| —A > l/(4r 2 ), which 
yields the general spectral bound 



E > min 

r>0 



p2 



We showed in Ref. [19( that general upper and lower estimates for all the discrete eigenvalues can be expressed in 
this form (for suitable P), with the aid of envelope theory [12, [H, [25|. If the exact eigenvalues of H are written 
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E(Z, (3, g), then these spectral functions are monotone in each parameter, decreasing in Z and q, and increasing in (3. 
Thus 

dE dE , dE 

dZ <0 > 5^ >0 ' and ^ <0 - 



Moreover, the following general scaling law is obeyed: 



1 




E(Z, 13, q) = Z 2 E(l, Z(3, q) = pE(Z(3, 1, q). (2) 

In the limit as q — > oo, the potential descends to the cut-off Coulomb potential Vqo given by 

-#, if r </3; 

lim V(r) =Foo(r) 

g — >oo 

This much is known for the whole class of problems. 

Mehta and Patil [4| have presented analytical solutions for the s-state eigenvalues corresponding to the V\ potential. 
Patil [5| has also discussed the analyticity of the scattering phase shifts for two particles interacting through the 
potentials V q with q = 1 and q = 2. Singh et al [1] have reported a large number of eigenvalues for the states Is to 
4/ corresponding to V\ and V2 for a fixed value of Z; these values were obtained by the numerical solution of Eq.(l) 
for Z = 1. The scaling law ([2]) extends their application to other values of Z. Exact bound-state solutions of Vi have 
been considered earlier 0, E3] , in which only a limited number of states with a specific choice of t = ... 3 have 
been treated. To our knowledge, no such study for q > 2 has been reported so far. 

The principal results of the present paper are as follows. We show that Schrodinger's equation for the two physically 
important cases, q — 1 and q — 2, can both be solved analytically, in terms of the Heun differential equation and the 
generalized Heun equation, respectively. Interestingly, in a suitable formulation, we show that exact solutions exist 
in which a factor in the wave function becomes a polynomial. It is shown that such solutions, and the corresponding 
exact eigenvalues, can be found explicitly by the use of an iterative method called the Asymptotic Iteration Method 
(AIM) |26l . I27I l28l . I29I l30l . l3ll |32[ HH, Hj, [35| . In cases where such polynomial solutions are not possible, AIM can 
also be used to construct approximate wave functions and eigenvalues. We summarize AIM in the next section, and 
we formulate the eigenvalue problem for general q > 1 in section 3. In sections 4 and 5 we discuss in detail the cases 
q = 1 and q = 2, respectively, providing tables of exact analytical solutions, where possible, and approximations in 
other cases. 



II. THE ASYMPTOTIC ITERATION METHOD 



The asymptotic iteration method (AIM) was first introduced as a technique capable of overcoming the computation 
complexity of solving Schrodinger's equation with singular potentials. It was first used as an approximation method 
for computing both eigenvalues and eigenstates, with the aid of a computer- algebra system. Because of its explicit 
iterative structure, the the method has subsequently been employed very effectively to obtain exact analytic solutions 
to many problems involving linear differential equations, especially those presented by non-relativistic and relativistic 
quantum mechanics. Recently, some progress has been made also in adapting the method to non-linear differential 
equations 36j. Perhaps because of its success in finding exact solutions, the strength of AIM as an approximation 
method has not been emphasized. For useful approximate solutions, the iteration sequence must be started with a 
suitable form for the wave function, and also a starting position r = r$. If these two seeds are well chosen, then 
either the sequence stops at an exact solution, or a non-ending iteration sequence is generated, which can then be 
terminated to yield a good approximation. Asymptotics usually provides a reliable guide to the initial form for the 
wave function; the starting point ro is then chosen so as to stabilize the resulting iteration process. The present work 
illustrates the effectiveness of AIM in finding both exact and approximate solutions of quantum-mechanical problems 
generated by soft-core Coulomb potentials. We shall now outline the method. Consider a differential equation of the 
form 

y" = X (r)y + s {r)y, (' = -7-) (3) 

ar 

where Ao(r) and So(r) are C°°— differentiable functions. A key feature of this method is to note the invariant structure 
of the right-hand side of (J3J) under further differentiation. Indeed, if we differentiate ([3]) with respect to r, we obtain 



y'" = X lV ' + siy 



(4) 



3 



where Ai = Aq + sq + Aq and s\ = s' + sqXq. If we find the second derivative of equation ([3]), we obtain 



2/ (4) = X 2 y' + s 2 y 



and 



respectively, where 



From (J6j) and (O we have 



y {n+1) = K-iy' + s„_iy 
y {n+2) = Ky' + s n y 



A„ = A(j_! + s n _i + A A„_i and s n = s' n _ x + s A„ 



An2/ (n+1) - A„_iy ( ™ +2) = S n y where <5„ = A n s n _i - A n _is n . 



(5) 



where A2 = X[ + Si + A Ai and s 2 = s[ + sq\i. Thus, for (n + l)" 1 and (n + 2)*' 1 derivative of n = 1, 2, . . ., we 
have 



(6) 
(7) 
(8) 
(9) 



Clearly, from J9J if y, the solution of ((3]), is a polynomial of degree n, then S n = 0. Further, if 5 n = 0, then <5 n / = 
for all n' > n. In an earlier paper [2^] we proved the principal theorem of the Asymptotic Iteration Method (AIM), 
namely 



Theorem 1: Given Ao and sq in C°°(a,b), the differential equation |3|) has the general solution 



y(r) = cxp I J a(t)dt 



C 2 + C! / exp / (A (r) + 2a(r))dr dt 



if for some n > 



s„-i 



A n A n _ 



(10) 



(11) 



For a given potential, the radial Schrodinger equation is converted to the form of equation ((3|). Then, so(r) and Ao(r) 
are determined and the s n (r) and X n (r) parameters are calculated from the recurrence relations given by equation 
(|8|). The energy eigenvalues are obtained from the roots of the termination condition of the method in equation 
or equivalently, 

d n (E;r)=X n (E;r)s n _ 1 (E;r)-X n . 1 (E;r)s n (E;r)=0, n=l,2,.... (12) 

where n represents the iteration number. For the exactly solvable problems, the energy eigenvalues are obtained 
immediately from the equation 

S n (E;r)=S n (E)=0, 

and very often the eigenvalue index is the same as the iteration number n. For problems that do not have exact 
solutions expressible in the form chosen, the approximations for the eigenvalues are taken to be the roots of this same 
equation for large iteration number. 



III. SOFT-CORE COULOMB POTENTIALS V q {r) 



In atomic units, the radial Schrodinger equation for the potential V q (r) reads 

1(1 + 1) Z 



1#_ 

2dr 2 



2r 2 



(r« + f3i)l 



tjj = Ei}). 



(13) 



4 



We may assume the solution of equation (fT"3|) . which vanishes at the origin and at infinity, is 

V>(r) = r I+lg-*(r«+^)« / (( r 9 + ^i) 



(14) 



Straightforward computation shows that f((r q +[3 q ) i ) is a solution of the second-order homogeneous linear differential 
equation 



r 2q-2( r q + /3 9 ) i~ 2 /"((r 9 + £«)«) + 



(g-l)^r 9 - 2 (r 9 + /3 9 )^ 



i-2 



+ 2vr q ~ 2 {r q + /3«)« -1 - 2kr 2q - 2 (r q + (3 q )«' 



/'((r« + /?*)*) 



% - l)/3 q r q - 2 (r q + f3 q )*- 2 - 2kvr q - 2 (r q + f3 q )^ 1 



+ k 2 r 2q - 2 (r q + P q )-«- 2 + 2Z{r q + f3 q y* + 2E 



f((r q + P q )*) = 0. (15) 



where, for simplicity, we write v = I + 1. From equation (|15j) . we obtain, after some simplification, 



f"(ir q + (3 q )~«) + 



(q - l)P q | 2^ 

r 9( r 9+/39)l r9(r9+/?9)9 

k(q - l)P q 2kv 



2k 



r q (r q + j3 q )« r q (r q +(3 q )~? 



+ A: 2 + 



2Z 



2£ 



/((r* + /3«)«) = 0. 
(16) 



If we denote x = (r q + /3 q ) * , we may now write (fTTj)) in a more compact form, as 
" (q ~ l)P q , 2v 



fix) 



xix q -P q ) x 1 - q (x q -P q ) 

k(q-l)j3 q 2kv 



2k 



fix) 



2Z 



2E 



X(x q -P q ) X 1 - q (x q ~P q ) X 3 - 2q (x q ~ q ) 2 ~« X 2 - 2q (x q ~ f3 q f~« 



fix) = 0, (17) 



where differentiation is now with respect to \. This equation is now in a form suitable for the application of AIM 
with: 



( \ M = (g- 1 )/ 3 " | 97, 



2fci/ 



x(x 9 -/3") " T " X 1_, (x 9 -/3 9 ) 



2Z 



2i' 



(18) 



In the next two sections we shall study the most important cases for physical applications, namely q = 1 and q = 2. 



IV. BOUND STATES OF SHIFTED COULOMB POTENTIAL Vi(r) = -Z/(r + 0) 

The problem of determining the energy eigenvalues and eigenstates of the shifted Coulomb potential 

Z 



Vi(r) 



13 > 



(19) 



has been of some interest in the past. As mentioned in the introduction, the potential may serve as an approximation to 
the potential due to a smeared charge distribution, rather than a point charge, and may be appropriate for describing 
mesonic atoms [37|. Using equation (Tj"4")) . the solution of Schrodinger equation with the shifted Coulomb potential 
V x {r) is 



iP(r)=r l+1 e- k xf( X ), X = r + p, 



(20) 
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where f(x) is the solution of differential equation 



fix) 



2v 



2k 



fix) 



2kv 



or, for k 2 = —2E, we have 



fix) 



2v 
X^P 



- 2k 



fix) 



2kv 
X 3 ^ 



2Z 

X 

2Z 

X 



2E 



fix) = 



fix) = o- 



(21) 



(22) 



By comparing this equation with equation l|18[). we can easily write the initial \o(r) and Sq(t) values, 

A (r)^-^ + 2fc, 
Sn ( r ) = Me. - M 

By means of AIM sequences, Eq. ([8]), we may calculate A n (r) and s„(r). Equation (fT2| then gives 
<?! = => (Z - -ku- k)x 2 - 2Z(i{Z -kv- k)x + Zb(Zf3 -k0-u)=O. 



(23) 



(24) 



From this equation, for k = Z/[y + 1), we have a condition on the parameter (3, namely Z(3 = v + 1. In general, we 
have in terms of the iteration number n (or the degree of the polynomial solution f(x) °f (|22j) l 



/c 



n = l,2,. 



(25) 



with conditions on the parameter /3 reported earlier [19(. For other values of k and (3, not necessary obeying these 
conditions, the eigenvalues are then computed by means of the termination condition (|12[) . namely S n (r; E) = 0. As 
we mentioned above, the computation of the eigenvalues by means of (|12|) should be independent of the choice of 
r. However, for certain values of ro, we may encounter oscillations of the computed roots and values that seem to 
diverge, presumably owing to rounding and computational errors in the algorithms used. In practice, this problem is 
avoided by re-choosing ro. In Table [H we show some results of these calculation obtained by means of the computer 
algebra system Maple. This environment allows numerical calculation of the roots of equation (|12[) with arbitrary 
chosen precision. In order to accelerate the computation we have written our own root-finding algorithm instead of 
using the default procedure Solve of Maple. In Table HJ we also report the iteration number along with the initial 
ro used to obtain the eigenvalues accurate to the number of decimal places recorded. For much smaller values of the 
parameter 0, a large number of iterations are usually needed, along with careful adjustment of ro- 



A. Analytic solutions for Vi(r) = —Z/(r + P) and the confluent Heun equation 



It is interesting that the polynomial solutions of f(x) of equation (f22|) are completely captured by AIM through the 
relations governed by k and (3 reported in Table 1. In other words, no other polynomial solutions are possible. In 
order to confirm this claim, we now express the solution of the differential equation (f2"2")l in terms of the confluent 
Heun function. From Appendix I, it is clear that the differential equation (j2"2"| is a special case of the confluent Heun 
differential equation (|69[) . Indeed with simple transformation z — xl Pi equation (|22[) can be written as 



2v 



dz 2 

Further, the transformation z = 1 — t yields 



-2k(3 



dt 2 



2v 
T 



2k(3 



df_ 

dz 



df_ 

dt 



2kv(i 2Z(3 



1 



2kv(3 2Z(3 



f = o. 



f = 0, 



t t-1 

which is clearly a special case of Exact solutions are given (see equations ([70)1 and (fTTj) ) by 

f{t) = He(2kp, 2v - 1, -1, -2Z/3, ~,t). 



(26) 



(27) 



(28) 
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TABLE I: Approximate solutions of the radial Schrodinger equation (I13|) for q = 1. 

(3 I Energy TV ro 

200 -0.003 653 168 9 14 210 

1 -0.003 169 532 8 17 210 

2 -0.002 798 561 8 17 210 

3 -0.002 498 271 8 17 210 

100 -0.006 742 076 7 17 120 

1 -0.005 633 737 9 19 120 

2 -0.004 812 493 4 23 140 
3 -0.004 168 784 4 13 150 

50 -0.012 194 692 6 16 75 

1 -0.009 717 588 4 18 100 

2 -0.007 962 796 3 13 100 
3 -0.006 643 881 8 16 110 

35 -0.016 388 672 4 20 65 

1 -0.012 685 924 3 16 65 

2 -0.010 135 716 7 12 75 
3 -0.008 268 680 7 14 90_ 

20 -0.025 669 937 8 38 65 

1 -0.018 846 220 6 30 65 

2 -0.014 387 207 9 22 65 

3 -0.011 278 293 3 22 75 
10 -0.043 438 719 3 56 45 

1 -0.029 446 515 7 49 55 

2 -0.021 024 301 6 43 65 

3 -0.015 576 600 1 42 75 



For the polynomial solutions of fit), the parameters of He(2kf3, 2v — 1, — 1, — 2Z/3, i,t) have to obey the conditions 
of equation (fT2)) . The first condition of ([72f immediately yields 

k - (29) 

where N is the degree of the polynomial solution. This is in complete agreement with the results of AIM as given 
by While the second condition of (|72")1 now yields the conditions under which these polynomial solutions are 

possible, namely, 

A N+1 (2kvp) = 0, (30) 

where Ajv+i(2fc^/3) is given by means of the tri-diagonal determinant of Table ILT1 Since t = 1 — z = 1 — $ = 1 — ^42 = 
— -g, the Taylor series expansion /jv(i) becomes 

Mt) = He (j*L t 2 V - 1, -1, -2Z0, \ - I) = £ u n (J*L, 2 V - 1, -1, -2Z0, i) (-l) " . (31) 
Therefore, we have 

/i(0 = i + T f : 7; (32) 

for 

A 2 (^£) ^3-^-1 = (33) 
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and 



if 



Further 



Zr Z{-v 2 + - 4)t/ + Z(3 - 4)r 2 
/2lrJ 2 + t/ + + 2) 2 (2//+l) 



A 3 = Z 2 (v + l)/? 2 - 3Z(i/ + 2)(iy + l)/3 + (2u + l)(v+ 2) 2 = 0. 



(34) 



(35) 



Zr Z(-v 2 + (ZP - 6)1/ + ZP - 9)r 2 
Mrj_ + 3 + ^ + /% + 3) 2 (2*/ + l) 

Z(2^ 4 + (19 - 3Z/?)z/ 3 + (Z 2 /? 2 - 21Z/3 + 63)i/ 2 + (3Z 2 /? 2 - 45Z/3 + 81)z/ + 2Z 2 /3 2 - 27Z/3 + 27)? 



if 



A 



/? 2 (6i/ + 3)(i/ + 3) 3 (j/ + l) 



4 J = Z 3 (v + 2)(i/ + I)/? 3 - 6Z 2 (i/ + 3)(i/ + 2)(i/ + l)/5 2 + Z(lli/ 2 + 28v + 15) (i/ + 3) 2 /3 



(36) 



3(i/ + l)(v + 3) 3 (2^ + 1) = 0, 



(37) 



and 



f( x = 1 . gr , Z(-^ 2 + (^-8)^ + Z/3-16)r 2 
/4lrj 4 + i/ /3(r/ + 4) 2 (2^+l) 

Z(2i/ 4 + (25 - 3Z/3V 3 + (Z 2 (i 2 - 27 Z [3 + 108)z/ 2 + {3Z 2 (3 2 - 72Zf3 + 17 6) v + 2Z 2 (3 2 - 48Z(3 + 64)£ 3 



Zr 4 



+ 6 /3 3 (, + 4) 4 (2,;i)(, + l)(2, + 3) ^ + - 105) ^ + ( - 6Z2/?2 + 163Z/? - ™»* + (918Z/? + 

- 2448 - 66Z 2 /3 2 )i/ 3 + (-252Z 2 /? 2 - 4128 + 2408Z/3 + 6Z 3 /3 3 )i/ 2 + (-3072 + 2848Z/3 - 384Z 2 /? 2 + 11Z 3 /3 3 )i/ 

- 768+ 1152Z/3 + 6Z 3 /? 3 - 192Z 2 /? 2 )) (38) 



if 



A 5 [j^jj = Z A {v + 3)(i/ + 2)(z/ + l)/3 4 - 10Z 3 (i/ + 4)(z/ + 3)(i/ + 2)(v + l)/3 3 + Z 2 (35^ 3 + 195^ 2 + 328j/ + 162) (i/ + 4) 2 /? 2 

- Z(50i/ 3 + 231^ 2 + 313^ + 126) (t/ + 4) 3 /3 + 6(2u + 3)(2n + l)(i/ + l)(i/ + 4) 4 = 0. (39) 
Similarly, for the higher polynomials, in general, 

A N+ x(2kv/3) = 

where /S.M+i{2kv0) is given by Table [TTJ 

All the N + 1 roots of this determinant are real and distinct [38| ■ It may be possible in future to obtain such results 
by studying the correspondence between the iteration details of AIM and the structure of the resulting tri-diagonal 
determinant. 



V. BOUND STATES FOR THE SOFT-CORE COULOMB POTENTIAL V 2 (r) = -Z/^/r 2 +(3 2 

For q — 2 and v = 0, we have using Eg. ([17]) 



f"(x) 



2 



2^X 

x(x 2 - P 2 ) x 2 - P 2 



2k 



fix) 



k/3 2 2kv\ 



2Z X , 2E X 2 



x(x 2 -P 2 ) x 2 ~P 2 x 2 -F x 2 -P 2 



fix) = o. 

(40) 
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TABLE II: The determinant A N +i(2Zj3). Here k = Z/(N + v) 

2kvf3 2u 

2Nk(3 2ku0 — (2i/ — 2fe/3) 2(2v + 1) ... 

2(N - l)k{3 2kv(3 - 2(2f + 1 - 2fc/3) ... 



2kv{3 — (TV — 2)(iV +2;/ - 2fc/3 - 3) (JV - 1)(JV -f 2f — 2) 

Akf3 2kv(3 - (N - 1)(JV + 2u - 2kf3 - 2) N {N + 1v — 1) 

2fc/3 2fci//3 - JV(JV + 2u + 2kf3 - 1) 



From this equation we have 



(i-/3W(x) + 



P 2 X ' 3 - 2k + 2k(3' z x 



2, ,-2 



/'(x) 



/(x)=0, (41) 



- k/3 2 X ' 3 + k 2 + 2ZX' 1 + 2E- k 2 f3 2 X ~ 2 

as given earlier by Liu & Clark |l7| for — 2E = k 2 . For the purpose of applying AIM, we may note for v = I + 1 and 

2E = -k 2 

that 



f"(x) 



2k 



2 



2v X 



X (x 2 -P 2 ) x 2 ~(3 2 
This equation can be further simplified, yielding 



fix) 



kf3 2 



2kvx 



k 2 [3 2 



2Z X 



fix) = 



2fc+ l_(2« + l)x 



X x 2 ~(3 2 



fix) 



xix 2 ~P 2 ) X 2 ~P 2 X 2 ~P 2 X 2 ~(3 2 
k {{2v + l)k-2Z)x + k 2 (3 2 ^ 



fix)- (42) 



X 



x 2 -P 2 



fix) 



fix) = 



2k+ i_ ("+V2) («+V2) 



X X-/3 X + 

Thus we may now apply AIM with 



fix) 



(u+ l/2)k- Z + k 2 /3/2 (y + \/2)k- Z -k 2 (3/2 k 



x-P 



X + P 



u x x-P x+0 ' 



so 



(v+l/2)k-Z+k 2 p/2 (v+l/2)k-Z-k 2 p/2 _ fe 
x-P + X' 



Using AIM sequence ([5]) we can thus produce a set of exact polynomial solutions given by 



k = 



v + n 



E n t — — — 



1 Z A 



2(n + v) 2 ' 



(43) 

fix)- 

(44) 
(45) 

(46) 



where n = 1, 2, 3, . . . is the iteration number used by AIM, along with the corresponding conditions on the potential 
parameter (3 reported in Table 11111 



The asymptotic iteration method can also be used to compute approximate eigenvalues of Schrodinger equation (|44|) 
for different explicit values of the potential parameters Z and (3. We use Ao and sq given by (|45p . and the numerical 
results, including the iteration number and suitable ro, are reported in Table pvp . It is important to note for the 
case of (3 = 200, the eigenvalues for I = 0, 1, 2, 3 are squeezed together in space less than 0.001, and this makes it hard 
for AIM to function correctly. In such cases, it is therefore recommended to use the scaling law as given by ©. For 
example, we know using the scaling law to show that 

£(1,200,2) = -^£(4,50,2). 

Thus instead of computing £(1,200,2) we may compute E(A, 50,2) then divide the result by 16. Indeed, for I = 2 
and I = 3, AIM yields for £(4,40,2), the values -0.062 403 136 6 and -0.057 981 633 1, respectively. From which 
we can now compute the eigenvalues E(l, 200, 2) as reported in Table (|IVp. 
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TABLE III: Polynomial conditions on the potential parameters for the existence of exact solutions of the radial Schrodinger 
equation (I44p for q — 2. 



Conditions on /3 



z 

u+l 



z 

v+2 



z 

v+3 



z 

I/+4 



z 

v+5 



z 

1/+6 



Z 
v+7 



Z 2 2 - 2{y + l) 3 = 

Z 4 /3 4 - 6Z 2 (v + 2) 3 /3 2 + 4(2u + + 2) 5 = 

Z 6 /3 6 - 12Z 4 (V + 3) 3 /3 4 + AZ 2 (llu + 18) (i/ + 3) 5 /3 2 - 24(2^ + l)(i> + l)(v + 3) 7 = 

Z 8 /3 8 - 20Z 6 (i/ + 4) 3 /3 6 + 20Z 4 (7^ + 19) (i/ + 4) 5 /3 4 - 8Z 2 (50u 2 + 193v + 182) + 4) 7 /3 2 + 96(2^ + 3)(2i/ + 1)0 + l)(i/ + 4) 9 = 

Z 10 /3 10 - 30Z 8 (k + 5) 3 /3 8 + 20Z 6 (17^ + 64) (i/ + 5) 5 /3 6 - 24Z 4 (75^ 2 + 456^ + 685)0 + 5) 7 /3 4 
+16Z 2 (274z/ 3 + 1821j/ 2 + 3942^ + 2765)0 + 5) 9 /3 2 - 960(2^ + 3)(2v + 1)0 + 2)0 + 1)0 + 5) 11 = 

^12^12 42 Z 10 (v + 6) 3 /? 10 + 140Z 8 (5f + 24)0 + 6) 5 ,3 8 - 168Z 6 (35i/ 2 + 288t< + 588)0 + 6) 7 /3 6 

+16Z 4 (1624i/ 3 + 16338^ 2 + 54099^ + 58860)0 + 6) 9 ,3 4 - 288Z 2 (196i/ 4 + 1956i/ 3 + 7157^ 2 + 11330^ + 6528)0 + 6) n ft 2 
+5760(21- + 5)(2i> + 3)(2i/ + 1)0 + 2){u + 1)0 + 6) 13 = 

Z 14 f3 14 - 56Z 12 (v + 7) 3 /3 12 + 56Z 10 (23t< + 134)0 + 7) 5 /3 10 - 112Z 8 (140^ 2 + 1447^ + 3719)0 + 7) 7 /3 8 

+16Z 6 (6769i/ 3 + 90363^ 2 + 398730i/ + 581140)0 + 7) 9 /3 6 - 64Z 4 (6566^ 4 + 96071^ 3 + 520284^ 2 + 1234091^ + 108054)0 + 7) xl /3 4 
+576Z 2 (1452^ 5 + 201121/* + 109039^ 3 + 288179^ 2 + 370287^- + 184667)0 + 7) 13 /? 2 
-80640(2^ + 5)(2v + 3)(2i/ + 1)0 + 3)0 + 2)0 + 1)0 + 7) 15 = 



TABLE IV: Approximate solutions of the radial Schrodinger equation (|13[) for q = 2. The values marked * are computed using 
a scaling law, as discussed in the text. 



P I Energy N ro 

200 -0.004 502 854 6 23 3 

1 -0.004 193 071 3 21 3 

2 -0.003 900 196 f 16 2 

3 -0.003 623 852 l 1 " 16 2_ 

100 -0.008 629 777 4 21 3 

1 -0.007 800 132 8 18 3 

2 -0.007 035 193 2 21 3 
3 -0.006 332 729 19 3_ 

50 -0.016 260 721 2 23 3 

1 -0.014 088 375 1 19 3 

2 -0.012 158 711 5 22 3 
3 -0.010 458 421 2 19 3_ 

35 -0.022 334 284 2 23 3 

1 -0.018 810 954 4 23 3 

2 -0.015 761 480 4 22 3 
3 -0.013 152 440 8 21 3_ 

20 -0.036 198 545 7 25 3 

1 -0.028 830 115 2 24 3 

2 -0.022 787 158 5 23 3 

3 -0.017 928 611 8 24 3 
10 -0.063 738 918 2 27 3 

1 -0.046 199 039 27 3 

2 -0.033 158 588 8 32 4 

3 -0.023 806 736 2 37 5 
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A. Analytic solutions for Vb(r) = —Zj \/r' 2 + f3 2 and the generalized Heun equation 



To analyze the analytic solutions of the differential equation (|44[> . we first let \ — j3z. We then have 

'{v + \/2)kj3- Z(3 + k 2 f3 2 /2 (v + l/2)fc/3 - Z(3 - k 2 [3 2 /2 kf3 



f" = 



z z — 1 z + 1 



Z - 1 



z + 1 



or, in the more appropriate form, 



2fc^- 



1 (w + 1/2) (w + 1/2) 



z - 1 



z+1 



f3k + f3 2 k 2 z + 2(3(kv - Z)z 2 
z(z-l)(z + l) 



/ 
(47) 



(48) 



This equation can be regarded as a special case of the generalized Heun differential equation. According to Schafkc 
and Schmidt [401 ], the generalized Heun equation has the form 



f" 



1 - Mo , 1 - Mi , 1 - M2 



z — 1 z — a 



fa + 0iz + (3 2 z 2 



z(z — l)(z — a) 



/ = 



(49) 



in which a £ C\{0, 1} and Mo, Mi> M2> <2> A), A) /?2 are arbitrary complex numbers. For a ^ 0, the generalized Heun 
equation (|49]) has three regular singular points, at 0,1, a with exponents (i.e. the roots of the indicial equations) 
{0, Mo}i {0,/Xi} and {0, m 2 }; respectively, in addition to one irregular singular point at infinity. These can be shown 
to be given by: 



1. For the singularity z = 0, we have 



lim z I a + 

z^Q 



1 - Mo , 1 - Mi , 1 ~ M2 

z z — 1 z — a 

/? + /3iz + /3 2 z 2 
z(z — l)(z — a) 



lim z 2 

2^0 



= 1 - Mo 
= 



and the indicial equation then reads 

s(s - 1) + (1 - fi )s = ^ s(s - Mo) = 

with simple roots {0, Mo}- 
2. For the singularity z = 1, we have 



lim(z — 1) [a + 



1 - Mo , 1 - Mi , 1 - M2 

z z — 1 z — a 

I3 + Piz + /3 2 z 2 



lim (z — 1) 
>i 



r(z — l)(z — a) 



= 1 - Mi 
= 



and the indicial equation then reads 

s(s - 1) + (1 - mi)s = s(s - mi) = 

with simple roots {0, Mi}- 
3. For the singularity z =a, we have 

1 - mo . 1 - Mi , 1 - M2 



lim (z — a) I a + 



z - 1 



z — a 



lim (z — a)" 

z— >a 



2 ^ /3 + /3iz + /3 2 z 2 
z(z — l)(z — a) 



1 - M2 




and the indicial equation then reads 

s(s - 1) + (1 - M2)s = s(s - M2) 

with simple roots {0, M2j-- 
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To find the exponent of the singular point at oo, we substitute z = — , Zx = -, dz\ — — \dz. Further 



and 



d 2 .f d df 



df df dz\ 
dz dz\ dz 

2 df \ dz\ 



__ = _ Z *V_ 

z 2 dz\ 1 dz\ 



1 df 



dz 2 dz dz dz\ \ 1 dz\ ) dz 
Substituting this in we obtain 

' 1 - Ml 



-2zx 



df 2 d 2 f \ , 2 \ „ i df 4 d 2 f 



Zl dzx + l dzl 



a + (1 — na)zx + zx 



l-zx 



Z\ 



1 - Mo 
1 — azi 



/3 z^ + / 3 1 z 1 +/3 2 
zi(l - zi)(l - azi) 



/ = 0. (50) 



Thus we have for a ^ that z = oo is an irregular singularity of Poincare rank 1, since the term z k+1 (a/zf) = a/ 
for k = 1. Clearly then for the Heun equation (j4"§)) with a = 0, the singular points 0,1, a and infinity are regular 
singular points. This is the main difference between the Heun differential equation (a = 0) and the generalized Heun 
equation as given by (|4*§)) . 

In order to find a power series solution of (|49[) about the singular point z — 0, i.e. valid for \z\ < 1, we notice first 
that equation (j4"5)) can be written as 



(z 3 ~{l + b)z 2 +az)f" 



az 3 + (3 - a(l + a) - Mo ~ Mi — ^2)z 2 + (-2 + /i + M2 + (-2 + Mo + Mi)a + aa)z 



+ (1 - /i )a 



/' + 



fa + faz + foz* 



/ = 0. 



(51) 



Using 



/(z) = ^ Cfe z fe , £ = 5> fe z*-\ ^XX*- 1 ^ 

fc=0 fe=l ' fc=2 



k-2 



and substitute in (|5ip . we obtain, after some simplification, 

oo oo oc oo 

^2(k - l)(fc - 2)c fe _iz fc - (1 + a) ^ k(k - l)c k z k + a^ k(k + l)c k+1 z k + a^(k - 2)c fc _ 2 ; 



fc=3 



fc=2 



fc=3 



+ (3 - a(l + a) - fiQ — fix - M2) y^(fe - l)c fc _iz fe + (-2 + Mo + /i 2 + (-2 + Mo + Mi)a + aa) kcfe* 



fc=2 



(1 -^o)a^(fc+ l)c fe+1 z fc +/5 ^c /£ z' £ +/3i^c fe _iz fe - 1 +S 2 ^c fe „ 2 z fc = 



(52) 



from which we obtain the four-term recurrence relation 



k=2 



[{k - l)(fc - 2) + (3 - a(l + a) - mo - Mi - Ma)(A; - 1) + Bx] c fc _i 
+ [-fc(fc - 1)(1 + a) + fc(-2 + ^ + M2 + (-2 + Mo + Mi)a+aa) + B ]c k 
+ [kk{k + 1) + &{k + 1)(1 - Mo)] c fc+ i + [a(k - 2) + B 2 ] c fc _ 2 = 0. 



(53) 



Replacing k by fc — 1, we obtain 



afc(/c - Mo)cfe 



(A; - 1) - (A; - 2)(1 + a) - 2 + mo + M2 + (-2 + Mo + Mi)a + cm + B 



Cfe-l 



(fc - 2) [ k - a(l + a) - Mo - Mi - M2 + -B 



a(k -3) + B 2 



Ck-3 = 



(54) 



with c_i = c_ 2 = 0. It should be noted that the Floquet solutions of Schafke and Schmidt 40] to the generalized 
Heun equation ([41?]) follows by using 



c k 



Tfc(Mo,Mi,M2,a, B ,Bx,B 2 ; a) 

r(fc + 1 - Mo)r(fc + 1) ' 



t_i = t_ 2 = 



12 



to obtain, using (|53p . the four-term recurrence relation 

(ft — 1) I ft — fxo — (i2 + (k — Mo - Mi)& — aa ) — _B 



ar fc 



Tfc-l 



+ (ft-2)^-ft + a(a+l) + /i +Mi+M2j - -Bi 
- [a(ft - 3) + B 2 ] (ft - 2 - /x )(ft - 1 - /xo)(fc - l)(fc - 2)r fc _ 3 . 
That is to say, in the notation of Schafke and Schmidt, 

r k =(j> l {k- l)rfc_i - 2 (fc - 2)r fe _ 2 + 03(fc ~ 3)T fc _3, 



(ft - l)(ft - 1 - /io)Tfc-2 



(55) 
(56) 



where 



0i (0 = £(£ + 1 - Mo - Mi) + + 1 - Mo - M2) - o£ - ^ 

<M0 = (£ + 1)(£ + 1 - Mo) + 2 - Mo - Ml - M2) - (1 + + ^ 

1 



MO 



a a 

-K +B 2 )(£ + 1 - Mo)(£ + 2 - Mo)(£ + + 2 )- 
a 

The condition for an TV-degree polynomial solution is equivalent to the equation cjy + 3 = in (|54[) . i.e. 

aN + B 2 = 0. 

If this condition is satisfied, the recurrence relation 
a(fc + l)(fc + 1 - Mo) c fc+i + 



(57) 



(58) 



(ft - 1)(1 + a) - 2 + /i + M2 + (-2 + Mo + Mi)a + ak) +B 



Ck 



(k - 1) ft + 1 - a(l + a) - /i - Mi ~ M2 + B 



Cfc-l 



a(fc-2)+B 2 



Cfc-2 = 



(59) 



for < k < N forms a system of N + 1 homogeneous linear relation in the N + 1 coefficients . This system has a 
non-trivial solution if and only if its determinant Ajv+i, given explicitly in Table V, is zero; that is to say, 



±N+1 



0. 



(60) 



Under the simultaneous fulfillment of the two conditions (|58|) and (|60"f , the solution of the generalized Heun equation 
P9"| reduces to a polynomial of degree N. In what follows, we use these conditions to compute the polynomial 
solutions of the differential equation . These results confirm our earlier findings using AIM and reported in Table 
(ED. 



Comparing Eq. (|48)l and Eq. (|49|) , we have 





-2, Ml 




1 


Mo 


= 2-f, M2 


~~ 2 


Bo 


= -/3ft, 


Bi = -f3 2 k 2 , 


T?2 


a 


= -2fc/3, 


a = -1. 





Condition (|58[) . immediately yields 



TV + ^ 



The polynomial solutions can be found by means of the condition (|6T)|) : for TV = 1, we have 



A 2 = 



-f3 2 k 2 (3k 
-2j3(kv-Z) 2v-P 2 k 2 



^k{(3 z k A - 2Z) = 0. 



(61) 



(62) 



(63) 



TABLE V: The determinant Ajv+i. 



B t -2 + MO + M2 + (-2 + MO + Ml + + B 2a(2-M0> OOO 

B 2 (3 - q(1 + a) - no - Ml - M2) + Bl 2( - 3 + MO + M2 + (-3 + MO + Ml + <*)&) + B 3a(3-Mo) 

B 2 + a 2(4 - a(l + a) - MO - Ml - M2) + «1 3(-4 + MO + M2 + (-4 + MO + Ml + <*)&) + B 4a ( 4 - Mo) 

B 2 +2q 3(5 - q(1 + a) - mo - Ml - M2) + b i 4( - 5 + Mo + M2 + (-5 + Mo + Ml + a)a) + B 
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From this, for k = Zj(y + 1), we have 



with polynomial solution given by 



Z 2 f3 2 - 2{v + iy = 



fl{z) = l+I^ z 



= 1 



v+l 

z 

v + l 

z 



1 + -Vr 2 +/3 2 - 

v+l 



For N — 2, we have 



-/3 2 fc 2 
-2(3(kv-Z) 




/?fc 
2v - /? 2 /c 2 3/c/3 
-2f3(kv - Z) - 2/3k 2 + 4z/-/3 2 fc 2 



from which, for A; = Z/(v + 2), we have 

Z 4 f3 4 - %Z 2 {v + 2 ff3 2 + 4(2i/ + + 2) 5 = 0. 

with polynomial solution given by 

2(3 2 Z 2 2 



h{z) = l + ~^Z r 



v + 2 I3 2 Z 2 -2{2v + l){2 + v) 2 ' 



= 1 + 



2Z 2 



i/ + 2 A P 2 Z 2 -2(2v + l){2 + v) 2 



X 



l + ^Jr^T0 2 + — 

v + 2 V (i 2 Z 2 -2{2v+l){2 + v) 2 



For N — 3, we have 



-(3 2 k 2 
-2f3(kv - Z) 





(3k 

2v - [3 2 k 2 3k/3 -3 

-2/3(kv - Z) - 2(3k 2 + Av-[3 2 k 2 5f3k 

-2/3(kv - Z) - 4/3k 6 + 6v-f3 2 k 2 



from which, for k = Z/(v + 3), we have 

Z 6 (3 6 - l2Z 4 {v + 3) 3 /3 4 + 4Z 2 (11^ + 18)(z/ + 3) 5 /3 2 - 24(2^ + l)(v + l)(v + 3) 7 = 0. 
with polynomial solution given by (note k — Zj{v + 3)) 



(64) 



(65) 



(66) 



(67) 



f 3 (z) = l + f3kz- 
= l + k X - 



2(3 2 k((3 2 k 2 - 6(1/ + l))(k(v + l)-Z) 



4kf3 3 (k(l + v)- Z){k(v + 2) - Z) 



(3 4 k 4 + 2fc(6fc - 5Z)f3 2 + I2(v + l)(2v + 1) f3 4 k 4 + 2fc(6fc - 5Z)f3 2 + I2(v + l)(2n + 1) ' 
2fc(/3 2 fc 2 -6(i/ + l))(fc(i/ + l)-Z) 2 4fc(fc(l + i/)-Z)(jfe(i/ + 2)-Z) 



:X 2 + ' 



f3 4 k 4 + 2k(6k - 5Z)f3 2 + I2{v + l)(2i/ + 1) A /3 4 fc 4 + 2fc(6/c - 5Z)/3 2 + I2{v + l)(2n + 1) 



2fc(/3 2 fc 2 - 6(1/ + l)){k{v + l)-Z) 
(3 4 k 4 + 2k(6k - 5Z)(3 2 + I2{y + l)(2v + 1) 
4fc(fc(l + v) - Z){k{v + 2)-Z) 



= 1 + k^Jr 2 + P 2 
+ 



P 4 k 4 + 2fc(6fc - 5Z)f3 2 + I2{v + l)(2n + 1) 
Similarly, higher degree polynomial solutions follow from the determinant Aat +1 = given by Table IV. 



(68) 
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VI. CONCLUSION 



The Schrodinger operator H = — ^A + V q (r) is of considerable interest in atomic physics. In this paper we have shown 
for the two most important cases q = 1 and q = 2, that the bound-state eigenequation Hip = Eip can be reduced to 
the problem of solving the confluent Hcun differential equation, or the generalized Heun differential equation. In the 
present context, the wave function can be written in a form in which analytical solutions with the correct asymptotics 
involve an unknown factor, which in certain cases can be shown to be a polynomial. A recently established Asymptotic 
Iteration Method (AIM) has been used to determine this factor, and the corresponding eigenvalue, explicitly in cases 
where the factor is a polynomial, and to construct accurate approximations in other cases. 



Acknowledgments 



Partial financial support of this work under Grant Nos. GP3438 and GP249507 from the Natural Sciences and 
Engineering Research Council of Canada is gratefully acknowledged by two of us (respectively RLH and NS). KDS 
acknowledges the Department of Science and Technology, New Delhi for the award of J.C.Bose National Fellowship. 
KDS and NS are grateful for the hospitality provided by the Department of Mathematics and Statistics of Concordia 
University, where part of this work was carried out. 



Appendix I: Confluent Heun differential equation 



In this appendix we summarize some facts about the confluent Heun differential equation, for more details and further 
properties, we refer the work of P. P. Fiziev [3^ | and the references therein. The confluent Heun's differential equation 
in simplest uniform shape is written as 



H" 



(3+1 | 7+1 
z z — 1 



H' 



z-1 



H = 0. 



(69) 



The equation has two regular singularities, at z = and z = 1, and an irregular singularity at z = oc. The constants 
fi and v are related to the constants a, (3, 7, 5, rj in the accepted notation of the confluent Heun function 



as follows 



H = He(a,(3,7,5,r),z) 



(70) 



(71) 



a(B+l) 

V = ~~2 - M - 



2 



The confluent Heun functions He(a, (3, 7, 6, rj, z) reduces to a polynomial of degree N > in the variable z if and only 
if the following two conditions are satisfied: 



S_ , £+2 

a "r" 2 



N + 1=0, 



(72) 



A N+1 (ji) = 

where Ajv+i(aO is a three-diagonal determinant given by: 



fi-qx 1(1 + /?) 
Na ii-q 2 +a 2(2 + (3) 
(N - l)a [i - q 3 + 2a 









M - qN-i + (N - 2)a {N - 1)(N - 1 + j3) 

2a fi-q N + (N-l)a N(N + (3) 

la jj. — qN+i + Na 



(73) 
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for 

q n = ( n -l)(n + /9 + 7) 

Under the simultaneous fulfillment of the two additional conditions (|72j) the confluent Heun function He{a, f3, 7, 5, r), z) 
reduces to a polynomial of degree N . 



[1] H.A. Bethe and E.E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms, (Dover Publications, New York 
1977). 

[2] M. Moshinsky, The Harmonic Oscillator in Modern Physics: From Atoms to Quarks, (Gordon and Breach, London, 1969). 
[3] K.D. Sen, V. I. Pupyshev and H.E. Montgomery Jr., Ad. Quantum Chem. 57, 25 (2009). 
[4] C.H. Mehta and S.H. Patil , Phys. Rev. A 17, 43 (1978). 
[5] S.H. Patil, Phys. Rev. A 24, 2913 (1981). 

[6] D. Singh, Y. P. Varshni and R. Dutt, Phys. Rev. A 32, 619 (1985). 

[7] H. De Meyer and G. Vanden Berghe, J. Phys. A: Math. Gen. 23, 1323 (1990). 

[8] A. Sinha and R. Roychoudhury, J. Phys. A: Math. Gen. 23, 3869 (1990). 

[9] F. M. Fernandez, J. Phys. A: Math. Gen. 24, 1351 (1991). 
[10] R.N. Chaudhuri and M. Mondal, Pramana-J.Phys., 39, 493 (1992). 
[11] M. Odeh and O. Mustafa, J. Phys. A: Math. Gen. 33,7013(2000. 
[12] O. Mustafa and M. Odeh, J. Phys. B: At. Mol. Opt. Phys. 32,3055 (1999). 
[13] C.A.S. Lima and L.C.M. Miranda, Phys. Rev. A 23, 3335 (1981). 
[14] J.H. Eberly, Q. Su and J. Javanainen, Phys. Rev. Lett. 62, 881 (1989). 
[15] Q. Su and J.H. Eberly, Phys. Rev. A 44, 5997 (1991). 

[16] M. Protopapas, C.H. Keitel and P.L. Knight, Rep. Prog. Phys. 60, 389 (1997). References therein. 
[17] C.W. Clark, J. Phys. B 30, 2517 (1997). 

[18] Y. I. Salamin, S. H. Hu, K. Z. Hatsagortsyan, C. H. Keitel, Phys. Rep. 427, 41 (2006). 
[19] R. L. Hall, N. Saad, K. D. Sen, and H. Ciftci, Phys. Rev. A 80, 032507 (2009). 

[20] M. Reed and B. Simon, Methods of modern mathematical physics II: Fourier analysis and self-adjointness, (Academic 

Press, New york, 1975). [The operator inequality is proved on p 169]. 
[21] S. J. Gustafson and I. M. Sigal, Mathematical concepts of quantum mechanics, (Springer, New York, 2006). [The operator 

inequality is proved for dimensions d > 3 on page 32.] 
[22] R. L. Hall, Phys. Rev. D 22,2062 (1980). 
[23] R. L. Hall, J. Math. Phys. 24, 324 (1983). 
[24] R. L. Hall, J. Math. Phys. 25, 2708 (1984). 
[25] R. L. Hall, J. Math. Phys. 34, 2779 (1993). 

[26] H. Ciftci, R.L. Hall and N. Saad, J. Phys. A: Math. Gen.36 11807-11816 (2003). 
[27] B. Champion, R. L. Hall, and N. Saad, Int. J. Mod. Phys. A 23 1405-1415 (2008). 
[28] M. F. Fernandez, J. Phys. A: Math. Gen. 37 6173 (2004). 

[29] H. Ciftci, R.L. Hall and N. Saad, J. Phys. A: Math. Gen.38 1147-1155 (2005). 
[30] H. Ciftci, R.L. Hall and N. Saad, Phys. Lett. A 340 388-396 (2005). 
[31] T. Barakat, Phys. Lett. A 344 411 (2005). 

[32] T. Barakat, K. Abodayeh and A. Mukheimer, J. Phys. A: Math. Gen. 38 1299 (2005). 

[33] T. Barakat, J. Phys. A: Math. Gen. 39 823-831 (2006). 

[34] O. Bayrak and I. Boztosun I, J. Phys. A: Math. Gen. 39 6955 (2006). 

[35] N. Saad, R.L. Hall and H. Ciftci, J. Phys. A: Math. Gen. 39 8477-8486 (2006). 

[36] N. Saad, R.L. Hall and H. Ciftci, J. Phys. A: Math. Gen. 40 10903-10914 (2007). 

[37] P. P. Ray and K. Mahata, J. Phys. A: Math. Gen. 22, 3161 (1989). 

[38] F. M. Arscott, Periodic Differential equations, (Oxford: Pergamon Press, 1964). The result is a consequence of Lemma 1, 
page 21. 

[39] P. P. Fiziev. larXiv:0904.0245 (2009). 

[40] R. Schafke and D. Schmidt, SIAM J. Math. Anal. 11 (1980) 848. 



